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1 Introduction 



Dispersionless integrable equations and hierarchies represent a particular 
class of integrable systems with a number of peculiar and remarkable prop- 
erties |l[]-[[n|. They arise in various problems of physics and mathematics 
from hydrodynamics and quantum field models to the theory of conformal 
mappings (see e.g. H,0-QIHf). 



Recently it was shown |fL9||-||21| that the dispersionless integrable hierar- 
chies are amenable to the quasi-classical 9-dressing method and that they 
are closely connected with the theory of quasiconformal mappings on the 



plane P2"|-[]2"l| . In these papers it was demonstrated that dispersionless 
integrable hierarchies are associated with the simple nonlinear <9-equation 

S- z = W(z,z,S z ). 

In the present paper we place this observation in a much wider setting. 
Namely, we consider a class of nonlinear problems on the plane which can be 
described by equations of the type 

S- z = W(z,z,S z ) + h(z,z), (1) 

where z,z£C and S, h and W are complex- valued functions. Such equations 
arise in several problems of hydrodynamics, electrostatics and quasiconformal 
mappings. We assume that nonlinearity W and source h are separated. We 
are looking for solutions of ([!]) in the form S = So + S, where So is determined 
by the source h. In this way the construction of solutions S of the problem 
([]]) is nothing but the dressing of the background solution Sq by the use of 
the quasi-classical <9-dressing method. 

Here we will concentrate on the study of properties of infinitesimal de- 
formations for the problem (p. We will show that these infinitesimal de- 
formations (J-p) obey universal (independent of the form of W) hierarchies 
of Hamilton- Jacobi (H-J) type equations, which in turn give rise to associ- 
ated hierarchies of dispersionless integrable systems. To derive these systems 
we use the quasi-classical <9-dressing method. Particular parametrizations 
of variations of sources h lead to different hierarchies of H-J equations and 
their associated dispersionless systems. Several concrete examples, like the 
dispersionless Kadomtsev-Petviashvili (KP) and the two-dimensional Toda 
lattice (2DTL) among others relevant examples, are considered. Equations 
arising in different gauges are also discussed. 
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2 General problem, parametrization of sources 
and the quasi-classical <5-dressing method 

So we will consider abstract nonlinear systems on the complex plane C de- 
scribed by an equation of the form 

S- z = W(z,z,S z ) + h(z,z). (2) 

The function W which defines the nonlinearity is assumed to be an analytic 
function of S z (i.e. independent of S z ). The inhomogeneous term h can be 
treated as an external source for the nonlinear system. 

Equations of the form (Q) arise in different fields of physics and mathe- 
matics. For example: 

1. Several problems of the plane motion of fluids [f25|l , generalized growth 
and Hele-Shaw problems |2^ . 

2. After differentiation with respect to z, equation @ reads 

S zg =(w(z,z,S z j) g + p(z,z), (3) 

where p = h z , so that it can be treated as the nonlinear Poisson equa- 
tion for the potential S. Such equation may arise in some special types 
of effective potential models in which W z and p describe an effective 
nonlinearity and a given external source, respectively. 

3. Under certain conditions solutions of (|2]) define quasiconformal map- 
pings of plane domains (see e.g.|27|]-[^]). Equations of the type (Q)arise 
also in the study of extremal problems for quasiconformal mappings 

(see 



Thus, the results obtained for equation (0) may have a wide range of 
applications. 

In our discussion we will consider the class of nonlinear problems (^) in 
which nonlinearity and sources are separated, i.e. we assume that there is a 
partition of the complex plane 

C = G Q \jG, G o f]G = 0, (4) 
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such that 



h(z, z) = 0, z E Go 

W(z,z,S g ) = 0, z E G. 
We will look for solutions of equation (§) of the form 

S = S + S, 



(5) 



where S is bounded on C and 



Sq.z — 0, z £ Gq 



S s = 0, zeG. 



(6) 



(7) 



The representation (0) imposes no constraint on the solution of 
In virtue of (@)-(0), equation (0) is equivalent to the system 



S 5 = W(z,z,So, z + S z ), zeC. 



(8) 



Thus, So is determined by the source h and can be considered as a back- 
ground (seed) solution of @ corresponding to W = 0. Hence, the whole 
procedure of construction of solutions of (§) (or the equivalent system (§D) 
is the dressing method, a well-known procedure in the theory of integrable 
equations (see e.g. |3Ofl-[0^|). A peculiar feature of the dressing method for 



(0)is that the corresponding <9-equation is a local nonlinear partial differen- 
tial equation. This (9-equation can be considered as the quasi-classical limit 
of the standard nonlocal 9-equation [I5|-[2I|.The method of construction of 



solutions within this quasiclassical <9-dressing method []19|-[pl| consists, ba- 
sically, in solving (H) or the quasi-linear equation for m = S z in the domain 



Go and the gluing with S in the domain G [ 20 1 . 

In this paper we will concentrate on the study of the properties of defor- 
mations of solutions of equation (0). Infinitesimal deformations (variations) 
S — > S + SS, generated by infinitesimal deformations h — > h + Sh of the 
source, are determined by the linear inhomogeneous Beltrami equation 



5S 



W'(z,z,S t }(6S) +5h, 



(9) 
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where W'(z,z,^) = 9W (*' Z & , Properties of the Beltrami equation are well- 
studied (see e.g.[03|,[p^|,|P3|). For our analysis we need two of them p3f : 



1. If fx, . . . , f n are solutions of the homogeneous Beltrami equation f z = 
fif z , then any arbitrary differentiable function F(f\, . . . , f n ) is a solu- 
tion too 

2. If the function \i satisfies \p\ < k < 1, then the only solution f z = fif z 
such that f s is locally LP for some p > 2, and such that / vanishes at 
some point of the extended plane C* is / = 0. 

These two properties provide a basis for the quasiclassical <9-dressing 
method p|, B The first 

one implies that, together with the infinitesi- 



mal deformations 5iS } (^S 1 , . . . , 5 n S, any differentiable function of them 

F(5\S, S2S, . . . , 5 n S) is also a solution of the Beltrami equation f z = W'f z 
in the domain Gq. Furthermore, under the conditions of the second property, 
if one is able to find a function F such that F(8\S, S2S, . . . , S n S) is bounded 
on C and vanishes at some point of C*, one has 

F(5 1 S,6 2 S,... ,5 n S) = 0. (10) 

The derivation of equations of the type (|10D is one of the main goals of the 
quasiclassical <9-dressing method. In this paper we aim to determine these 
kind of relations for the infinitesimal deformations of the problem (0). To 
do that one has to parametrize the source h and its deformations in one or 
another way. 

Having in mind the relation S 0:Z2 — h z — p(z,z), it is quite natural to 
choose p in the form 

N 00 

p = h z = J2Y.^ {n '°\z-z a ), (11) 

a=l n=0 

where 7 an and arbitrary complex constants (z a E G) and <5 <n '°' l (£) are 

the ^-derivatives of Dirac delta functions. Since 

Sq(z,z) = J J dz' A dz' la(z — z') p(z' , z') , 

one has 

N 00 

S (Z, Z)=J2 (toa -Z«)+J2 T J T L ^) ' ^ 

n=l \ Z Za > 



a=l 
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where t 0a = j a0 and t na = (-l) n (n - l)\j na , n > 1. 

Of course, all pole-type singularities in (|1^) can be obtained by coalescing 
logarithmic terms, so that only the S^' ^ terms in ( |TTD and the to a ln(z — z a ) 
terms in (|T2|) are, in fact, of fundamental importance. But it is convenient to 
add other terms to the source h from the very beginning instead of performing 
coalescence at the end. 

By considering an infinite number of points z a and by passing from the 
sum in ([11]) to an integral, one gets a source p of generic form. However, 
the simpler forms (|i~l~l) and (|T^) are much more convenient for performing 
calculations. The symbols t na and z a (a — 1, . . . , N; n > 0) are free pa- 
rameters, so that within the class of sources given by ( |TT1) the deformations 
(variations) of the sources are generated by variations of t na and z a . In this 
case 5h = and 5S = e|^, where e is an infinitesimal parameter and r is 
any of the parameters t na and z a . Hence, we have 

(S T )- z = W'-(S T ) z + h T . (13) 

In virtue of (Q) and ( p"2D one has 

St 0a = Hz - Za) + S t0a , 



(14) 



t>na Z — Z a 



and 



z- z a ^ [z - z a ) n+1 

n>l 

These expressions have singularities, but all of them are located in the domain 
G. 

Now let us consider a differentiable function of the type F(S tna , S Za , . . . ). 
It is a regular function on the domain Gq. Hence due to the first property of 
the Beltrami equation, one has 

F- Z = W'F Z , zeG . (16) 

In general this function has singularities in the domain G due to the singular 
terms of Sq t . But if we manage to cancel these singularities by a right choice 
of F, then we have 

Fg — 0, z e G, (17) 
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so that 



F- z = W'F Z , zeC. (18) 

In this case, as we are assuming that S is bounded in C, if F(S tna , S Za , . . . ) 
vanish at some point of the extended complex plane then from the second 
basic property of Beltrami equation we conclude 

F(S tna ,S Za ,...) = 0. (19) 

Equations like ( |T9| ) are our desired equations relating infinitesimal symme- 
tries of the system The next sections of the present paper are devoted 
to the derivation of hierarchies of equations of this type. As we shall see the 



form of the equations ( |19"D does not depend on the choice of W, so that they 
have a universal character. 

The problem (|2]) and the equations (|T9|) have an additional important 
property; namely, they are invariant under the gauge transformations 

S(z, z, t) -> S'(z, z, t) = S(z, z, t) + g(r), 

where g{r) is an arbitrary function. For concrete equations this property has 
been observed earlier in []I0| . This invariance implies that we can arbitrarily 
prescribe the value of 5* at any given point zq (except at one of the z a ), i.e. we 
can always normalize S by the condition S(z , r) = go(t na , z a ), in particular 
S(z ) = 0. At the point z = z a , one can choose the gauge S(z a ) = 0. As we 
shall see different gauges produce different, but related, equations of the type 
([T9|). It is also possible to give a gauge invariant formulation of equations 

&■ 

Before proceeding to the construction of concrete equations we would like 
to note the following: 

Remark 2.1 Equation (3) is the Euler-Lagrange equation with the La- 
grangian 



L = ^S z S g -Q(z,z,S z )+pS 

where ^ffiy^ = W(z, z, S z ). Equation (3) can be also represented in the 
form of the conservation law 
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{\s z )z + {\s- z - W(z,z, S z ) - h) z = 0. 

Formally this conservation law is connected with the gauge invariance men- 
tioned above. 

Remark 2.2 It is well-known that for the plane potential flows the ki- 
netic energy of a fluid coincides with the Dirichlet integral (see e.g. [34]). In 
our case the total "kinetic energy" ~ J f c S z Szdz A dz diverges due to sin- 
gularity of So. A regularized "kinetic energy" can be defined as the complex 
Dirichlet type integral 

D = \ [ f (OS A OS — 3S A 0S ) = \-. \ f(dSA*dS- dS A *dS ). 
2 J Jc 4« J J c 

Remark 2.3 Within the fluid mechanical interpretation of the prob- 
lem (0) the 7o a S °\z — z a ) terms in the source p or the logarithmic terms 
t 0a ln(z — z a ) describe sources or sinks of the fluid, while other terms describe 
vortices located at the points z a (see e.g. [0). In the electrostatic applica- 
tions the ~fo a S^°'°\z — z a ) terms or t Qa \n(z — z a ) terms, obviously correspond 
to point charged particles with charges t 0a . Other terms describe contribu- 
tions from electrical multipoles. Thus, the deformations we are discussing 
here are generated by the infinitesimal variations of the strengths and po- 
sitions of sources or sinks of fluid, charges of the point particles and their 
positions, and strengths and positions of vortices and multipoles. 

Remark 2.4 Equations of the type ( |i~9|) can be treated as the dispersion- 
less limit of certain dispersive integrable hierarchies. Under this limit wave 
functions ift of the fully dispersive case become ip = exp(^) where e — ► 0. 
Within such a connection the variables t na represent slow variables associ- 
ated with the standard flow variables x na = —. The variables t 0a can be 
considered as a quasi-classical limit of the Miwa variables. Indeed, in terms 
of Miwa variables, the undressed dispersive wave function ip has the form 
iPq{z) = ria=i(~ — 1) P " wri ere p a are integers (see p5| ). By considering large 
p a and by introducing t 0a via p a = t 0a /e (e — > 0), one gets 

; Ell^OaM^-^a) /Son 

ip ~ exp ~ expf — 1. 

e e 
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3 One point-like source case 

We begin by the simplest case of source p concentrated on a single point, i.e. 

00 t 

So = toHz-z a )+J2(7Z TK ( 20 ) 

n=l \ Z Z °-> 

Therefore 



(21) 



as _ 1 1 as 

dt n (z-z a ) n ' dtn ' 

From these equations we deduce 



and 



as 1 _as_ 
e 9t = e 8t ) ^ e c, (22) 



e-^W = i_ile-^)^, zeC, Zl eG, (23) 

2r 



where P(^) stands for the quasiclassical vertex operator 



1 rt 

-(;-,„)«_. (24) 



n=l 



Let us first consider variations of the the t n variables (time variables). By 
counting singularities of the derivatives Jp, one concludes that the simplest 
of the equations of the type fll4[) is of the form 

S t2 - (S tl ) 2 - vS h -u = 0, 

(25) 

v = -2S tl (z a ), u = S t2 (z ) - (S tl (z )) 2 - vS tl (z ), 
where zq is an arbitrary point zq 7^ z a . Other analogous equations look like 

n-l 

Si.-CSO"- j>(Sfc)* = 0, n ^ 2 ' ( 26 ) 

fc>0 
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where the functions Vk = Vk{t) are appropriate combinations of derivatives 
of S(z a ) and S(zq) with respect to t x . 

Equations (p6|) form a complete set of equations for infinitesimal varia- 
tions S tn of S, generated by variations of the times t n . They are equations 
of Hamilton- Jacobi type. 



An important property of equations (26) is that they have the same form 
for all functions W(z, z, S z ) describing nonlinearities in the problems (0). 

Remark 3.1 Within the quasi-classical <9-dressing method, the functions W 



play a role of quasi-classical <9-data fl9f,[pl]1. The independence of the form 



of integrable equations on the <9-data is a general feature of the dressing 
method. 

By construction, all the equations fl2~6f) are compatible and solvable by the 
use of the <9-problem (|2|) and, in turn, they determine a hierarchy of equations 
for the functions u, v and so on. The form of these equations depends on the 
gauge choice. 

There are two natural gauges. The first one is S(zo) = 0, (zq ^ z a ), in 
which u = and the hierarchy (p6|) becomes the Hamilton- Jacobi system of 
equations for the dispersionless modified KP (dmKP) hierarchy. The lowest 
member of which is IB6H 



v t3 + \v 2 v tl - \v tl d^{v t2 ) - \d^\v t2t2 ) = 0. (27) 



The second natural gauge is S(z a ) = 0. It leads to v = 0, and then (|26|) 
reduce to the system of Hamilton- Jacobi equations for the well-known dKP 
hierarchy (see e.g.p|-fl^fl). The dKP equation itself is 

-uu tl + 

Other gauges can be of interest too. For example, if we impose S(z ) = 
i.e. S(zq) = 0, we get a hierarchy of 2 + 1-dimensional Gardner (mixed 



ut 3 = ~ uu h + -d tl \u t2t2 ). (28) 



ZQ-Z a 

KP-mKP) equations. We shall refer to the hierarchy of equations ( P6|) and 
associated hierarchy for u and v as the d(KP-mKP) hierarchy. 

One can formulate equations in gauge invariant form. In such a formu- 
lation one has equations fl2"E| ) for the gauge invariant object S* := S(z) — 
S(z ), while equations for potentials Vk are formulated in terms of the gauge- 
invariant 



w := u + -d t v t2 - -v 
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This formulation provides us also with the dispersionless Miura transforma- 



tion 36 



Vl 1 2 



By proceeding along the same lines as in ||37|| , one can derive several 
important generating equations and addition formulae associated with equa- 
tions (|26|) . In virtue of (p3|) , both p := |j| and exp(— V(z\)S{z)) have a 
simple pole at z = z a , so that one gets 

p(z) - p(z ) = —U-vteW*) - e -^)S(zo)\ ; (29) 

Z\ — z a V / 

where z G C is an arbitrary point z ^ z a . 

Equation (|2§p is of fundamental importance for describing infinitesimal 
deformations. It is the generating equation for the hierarchy of equations 
(|2Bj). Indeed, by expanding both sides of ( |2UD in series on z x , one gets a 
system which is equivalent to (pp. Furthermore, by evaluating both sides of 
(|29|) at 2 = 2i, one gets 



p(z) - p{z x ) = — exp ( - V{z x ){S{z) - S(z a )j) , z G C, z x G G. 

Z\ — z a V / 

(30) 



Thus, by setting z = z 2 G G and using (p3|) and the skew-symmetry of the 
left-hand side of ([30]) under the interchange z\ ^ z 2 , one concludes that 



V(z x )(S(z 2 ) - 5(z a )) = V(z 2 )(S(z x ) - S(z a )). (31) 
Equation ( pi]) implies that there exists a function F such that 

S(z) = S{z a ) - D(z)F. (32) 

So one has 

p(z 2 ) - p{z x ) = e ^(^, ^ Z2 G & (33) 

(£q — Z\)\Z 2 — Z Q J 

Further, by summing up the relation (B3p for pairs of points (z 2 , z x ), (z x , z 3 ) 
and (z 3 ,z 2 ), one gets 

(zi - z 2 )(z 3 - z a )e v ^ v ^ F + (z 3 - Zl )(z 2 - z a )e v ^ v ^ F 

(34) 

+ (z 2 - z 3 )(^i - ^)e^^^ F = 0, z x , z 2 , z 3 G G. 
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It is the addition formula for the dKP hierarchy in the case when the sin- 
gularity is located at a finite point z a . One can get the usual form of the 
dKP addition formula (see e.g.|Il|,P2|,P!B|) by sending z a — > oo. Equation 



(|34|), which is gauge invariant, means that F = In TdKP-, where TdKP stands 
for the tau-function of the dKP hierarchy f/jHhflTll, g2|. 
From (E9|) one finds also that 



\n\p{z) - p(z )] - ln(z a - Zi) - V(zi)S{z a ) 



(35) 



In 



e -V( Zi )S(z) _ e -V(zi)S(z ) 



Zi eG,(i = 1,2,3). 



Acting on this equation by the operators T>(zk) and "D{z\) with different 
values of the indexes k, I, i, and summing up the corresponding equations, 
one gets 



^ e ik i{V(zi) - V(z k )) In 

i,k,l 



e -V{ Zl )S(z) _ e -V( Zl )S(z ) 



0, z 6 C, Zi e G, 
(36) 



where e^j is the totally antisymmetric tensor (eu3 = 1)- Equation (|36| ) is the 
generating equation for the hierarchy of equations for S(z, z, t). Indeed, the 
taylor expansion of its left-hand side leads to an infinite family of equations 
for S only. For example, in the gauge S(z a ) = 0, the first nontrivial order 
gives 



-5*1*3 — 4*5*2*2 + 2 



-5*2 _ 



(37) 



that is the S'-equation for the dKP equation pi| . In the gauge S(z ) = 
(zq ^ z a ), one has the corresponding equation for the dmKP hierarchy. Equa- 
tion (^) encodes all relations between infinitesimal deformations generated 
by variations of times t n . 

Now let us include also into consideration the variations of strength of 



the logarithmic term in (p0|). Using (|2l|) and (P2|), one readily gets 



S h = A u e- St ° + A 



10, 



S t2 = A 2 n e- 2St o + A 21 e~ St o + A 20 , 



(38) 
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where 



(39) 



A u = e 5t o(^), A 20 = S t2 (z ) - A\ x e- 2St ^ - A 21 e~ St ^\ 
A 21 = 2A 11 (%^) , A 10 = S tl (z )-A u e- s *oM. 

V / z=z a 

In general, 

n 

S tn = Y,A nk e- kSt °, (40) 



k=0 



where A n k(to,t) are certain functions. 

The hierarchy of Hamilton- Jacobi equations (HDf) is a basic one for the 



one-point-like source fl20|). Equations (|40"D imply a hierarchy of integrable 
equations for the coefficients A n ^. The lowest member of this hierarchy arises 
from the system of equations ( |38D and in the gauge S(zo) = looks like 

K+K~ eK =0, 

(41) 

t2 + 2^ 1 + (0 tl ) 2 + (^ o ) 2 = O, 

where (ft := S(z a ), tp := S z (z a ). 

Equation fl28|) allows us to derive an analog of the relationships (|29f)-(|34D 
for the logarithmic times. Indeed, by substituting ( p8[) into fl29"|) one obtains 

e -VS{z) _ e -VS(zo) + 1 e (P+P(*i))gQe) . 
^1 — 

(42) 



e -V( Zl )S(z) _ e -V(z 1 )S(z ) 



where we denote V = 



Similar to the case (|2T]|), equation is the generating equation for the 
hierarchy of equations (f40|). They are obtained by the expansion of the left- 
hand side of ( fED in Taylor series in Z\. Then, by evaluating ( pE2"D at 2: = 

one gets 

g-X>(5(*)-S(*»)) _ e -V(S( Zl )-S(z a )) + 1 c -2?( Zl )(S( 2 )-S( Za )) = Q) ^ 43 ) 

Zl — £ a 
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where 2 G C,Zi G G, which implies fl52|). Evaluating (^) at z = z 2 G G and 
using (|23|) and (|3^), one readily obtains the addition formulas for F: 

1 _ e VV(z 2 )F _ 1 ^Vjz^F _j Z 2 - Zl c D(z 1 )V(z 2 )F = q 



Z 2 — Z a Z\ — z a (z\ — z a ){z 2 — z a ) 

(44) 

where Zi,z 2 G G. Formulae (^0|) ,(|4lD-(^4|) incorporate those of the d(KP- 
niKP) hierarchy. It is easy to see that (fP|) implies the dKP addition formula 
(|34]) . It is straightforward to check that the function S, which obeys fllPD , 



solves also equations (|26| ) and, in particular, equations ( |38| ) imply (|25|) . Thus, 
it means that the logarithmic source is fundamental for the description of 
deformations of the problem (fj). 

Finally we consider variations of the positions z a where the source p is 
concentrated. For the pure logarithmic source S = t Q \n(z — z a ) and 

(45, 

Using (|22l) , one deduces that 

dS{z 



dz n 



+ Ae~ St ^ z) +B = 0, (46) 



where A = t expS to (z a ) and B = ~-§^(z )-Aexp (S to (z )), (z 6C,z / 
z a ). 

In the generic case ( |20| ) 

dS_ = t + nt n + 9S_ ^ 

dz a z z a \z z a ) n ^~^ dz a 

n>l 



In virtue of (ETT) one has 



AS = j°_ + y + Ag, (48) 

where 



2 Z n [Z z n 

n>l 



A3 ^£ + £ frt ._<L. (49) 

n>± 
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Using P7| ) and (f£8|), one concludes that 



<9S* 

— -AS + E = 0, (50) 

uz a 

where £ = + AS(z ), eC,z„/ z a . 

In particular, in the gauge S(zo) = 0, we have E = and 

This relation shows a connection of the variations of z a with non isospectral 
deformations for the dKP hierarchy. 

Remark 3.2 Equations (34), (36) and the dKP hierarchy are invariant 
under the scale transformations 





-4 = 


At„, n=l,2,..., 


— oo < A < oo 


F(t) - 


-> F(t' 


) = A 2 F(t), 




S(t) - 


■* 5 '(4 


) = AS(t), 




u(t) - 


- u'(0 


= u(t). 





Variation of the function F due to the infinitesimal variation of the times 
t n ( 5t n = et n ) is equal to 5F = e Y, n >i f n§^- 

Remark 3.3 Let us introduce the complex Dirichlet type integral 



2wi 




D G = — I lA dS A dS - dS o A 9S ). 



G 



For the d(KP-mKP) hierarchy one gets, using (20) with t = 0, the 
expansion S = J2 n>0 (z — z a ) n S n as z — > z a and the formula (32): 

Dq = — 2_j n ^nS n = ^ n 7fj~' 
n>l n>l n 

Thus, SF = sDq. An interrelation between the Dirichlet type integrals 
and r-functions for dispersionless hierarchies will be discussed in a separate 
paper. 
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4 Two-points source. Dispersionless Laplace 
and 2DTL hierarchies 



The case of a source p{h) concentrated in two different points z a and zp is 
rather rich. The domain G can be a disconnected set (G = V a {jT>p), where 
V a and Dp are two disks centered at z a and Zp respectively. Several types of 
different reduced hierarchies can arise. 
The function Sq has the form 

oo oo 

S = t 0a Hz -z a )+Y] Xn + t p \n(z -zp) + Y] - — . (52) 

^ {Z - ZaY ti(z~ ZpY 

Let us first consider variations of the variables x n and y n . Similar to the 
one-point case ( [26] ) one obtains two hierarchies of equations 

n 

S Xn =J2Uk(x,y)(S Xl ) k , (53) 

fc=0 



fc=0 

where Uk and Vk are certain functions. Each of these hierarchies gives rise to 
the dKP-niKP hierarchy in the variables x n and y n , respectively. Intercon- 
nection between both hierarchies is provided by the equation 

S Xl S yi + aS Xl + bS yi + c = 0, (55) 

where 

a = -Syi( z a), b = -S Xl (zp), (56) 

and 

c = -S Xl (z )S yi (z ) - aS Xl (z ) - bS yi (z ). (57) 

Here z eC, z / z a , Zp. 

Equation (^) can be treated as the quasi-classical limit of the well-known 
Laplace equation 

if) to + Aip^ + Btp v + aj) = 0. (58) 
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Indeed, by introducing slow variables X\ := e£, y\ := en, considering the qua- 
siclassical wave-function limit ip = exp(S/e) and proceeding to the limit e — > 
0, provided that A(f) -> a(xi,3/i), 5(f) -> 6(xi,3/i), C(f) -> c(x 1>3/l ), 
one converts ( p8| ) into fl55|) . 

So one can refer to the hierarchy described by equations d53|)- (^5|) as the 
dispersionless Laplace hierarchy It contains different interesting particular 
cases. In the gauge S(zq) = 0, the lowest member of this hierarchy is given 
by the equations 



S X iS yi + aS Xl + bS yi + c — 0, 



S X2 — {^Xi) ~t~ UiS Xl , Sy 2 — (Syi) "I" VlSy 1 , 



zp — z a 



(59) 



where U\ = —2S Xl (z a ), V\ = —2S yi (zp). Evaluating equations (j59|) at z 
z a , zp, one gets the following equations 

U X2 = —?—(U xl - V X1 ) + (U X1 ) 2 - 2U X1 V X1 , 



V X2 + {V Xi y = 24> xv (60) 

u Xl v yi - — L- v yi + —^—u xi + <j> yi = 0, 

z a — Zp Z a — Zp 



V, 



8/2 
2 



Uy 2 + (U yi J 



(z a - zp) y " :n 
2^ 



(V yi -U yi ) + (V m ) 2 -2V yi U yi , 



HI 



where U = S(zp), V = S(z a ), = § \ z=Za , = § \ Z=Z0 . 

Equations ( [59"D and (BO) are the lowest members of the dispersionless 
Laplace hierarchy. In terms of the variable t := x 2 + JJ2, equations (|60|) are 
equivalent to the dispersionless limit of the system of equations considered 

m 



m 



(equation (9.6)) and in |39 



39 



arc 



The equations discussed in 
the equations for the Davey-Stewartson wave function. One more case of 
interesting application of the dispersionless Laplace transform corresponds 
to the constraints S(z a ) = S(zp) = const. This constraint is compatible with 
equations (|53|) and (0) for odd n. The lowest member of such a hierarchy 
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is given by the equations 

S xl Sy 1 +(7 = 0, 

S X3 = (S X1 ) 3 + U 1 S X1 , (61) 
^2/3 = (Syi) 3 + ViS Xl , 

where U\ = 3d~^q Xl , V\ = 3d~^q yi . The corresponding equations for the 
potential q are 

q X3 = 3(q(dy 1 1 q xi )) xi , q m = 3(q(d-^q yi )) yi . (62) 
In terms of d t = d X3 + d y3 one has 

q t = 3{q{d~ 1 q Xl )) Xl + 3{q{d Xl 1 q yi )) yi . (63) 



Equation (^3|) is the dispersionless limit of the Nizhnik-Veselov-Novikov 
equation Pl-M 



qt = q Xl x lXl + q yiyiyi + ^>{.q{.d y ^q Xl ) Xl + 3(q(d^q yi 



yi ■ 



Notice that the dispersionless limit of the Veselov-Novikov equation was dis- 
cussed for the first time in f4|. 

Equations ( |6lD have an interesting application in physics. Namely, the 
first equation is, in fact, the eiconal equation from the geometric optics. 
Indeed, if we denote z — x%, z — yi, z = x + iy, it reads 

S 2 X + S 2 y = -Aq. 

Remark 4.1 For the multipoint sources case, the gauge S{zq) = const., z 7^ 
z a , zp is a symmetric and natural one. 
Remark 4.2 Equation fl55|) implies that 

s » - "trf (64) 

This relation shows that, using only one source z a , one can generate in- 
finitesimal deformations of S which have poles at other points (in this case 
at z = Zp). This fact, which was already discussed in Section 2, means that, 
in particular, it is convenient to introduce the corresponding times from the 
very beginning. 
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Let us include now the logarithmic times into consideration. Analogously 
to (f|(3), we have in the gauge S(z ) = 

n 

s Xn = j2 A * e ~ kStoa > 

(65) 



k=0 



where and Bj. are functions depending on the whole set of times. In 
addition there is an equation relating the derivatives of S with respect to the 
times t 0a and to/3 



e - s t Oa St o0 + Ae -s t0a + B e ~St og + c = Q } (66) 



where 



2/3 — z a 



B = — - — e ~ St ^\ C = -l-A-B. 



(67) 



Equations (p5|) -(|6q) describe the general two-points hierarchy of the in- 
finitesimal deformations. It contains a well-known dispersionless integrable 
hierarchy. Indeed, let us introduce new variables T and t such that t 0a = 
T + t, t 0a = T — t. Then one has 

e ±s t = (z_^A ±l e ±~s t _ (6g) 

\Z — Zp/ 

Using (|68|) , one obtains 

n n 

S Xn =J2Uke- kS \ S yn = Y,Vke kS \ n > 1, (69) 

k=0 fc=0 

for certain functions Uk and V*.. The lowest member of the hierarchy fl69| ) is 

S X1 = Ue- St - U, S yi = Ve St - V, (70) 
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where 



U 



1 



,s t (z a ) y 



-St(*/i) 



Z a ~ Zp Zp — z a 

Evaluating the left-hand sides of the equations flTOD at z = zp and z a , one 
gets 



XI 



£/3 — ^Q! 



1-e' 



where = S(zp), ip = S(z a ), or 

+ [7^ = 0, V X1 - VU t = 0. 

The system ( |7T| ) implies 



where 



= ln 



@xiyi + ^ e 



(z a - Z^) 5 



0. 



lg(C/F). 



(71) 

(72) 
(73) 
(74) 



Equation ( [73| ) is the well-known dispersionless 2DTL equation (see JTTj]) . Its 
associated system of Hamilton- Jacobi equations is provided by ( |69D [pH,|2l|. 
Note that equations flTUD imply the following equation for 5 



s xiyi + s xi s yi (^ e 5 t zryX ~~ ^ - 



(75) 



Now we will derive generating and addition formulae for the dispersionless 
2DTL hierarchy in the gauge S(zo) = 0. To this end we introduce the 
operators 



V a (z) :=£_, 

n=l 



Z - Z D 



dx n 



n=l 



A. 

dy n 



(76) 



One has 



e -23«(«i)S(z) = £ 5. e -»«(a)S(2) 



(77) 



z - Zp 
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where zeC, Zi,z 3 G G. Using ( [77|) one obtains for p a := and p/3 := J^- 

Pa(z)-Pa(z 1 ) + — e -".W(SW-SW) = 0, 2 6 C, 2i 6 G, 

— z a 

(78) 

pp{z) - pp{z 2 ) + — - — e -%(^)(s(^)-5(^)) = o, z g C, z 2 e G. 

z 2 — 2/3 



By proceeding as in the one-point case, one finds that equations (|78|) 
imply that 

S(z) = S(z a ) - V a (z)F a , S(z) = S{z p ) - Vp{z)F p , z G G, (79) 

where F a and Fp are functions depending on the times only, which satisfy 
addition formulae of the form fl34|). 

On the other hand, using (^) one gets the identities 

Pa(z)-p a (z p ) = - J_ e -^W-^)), 

2 Q — 2/3 



(80) 



i p -D{S(z)-S{z fj )) 



P(3{Z) -Pp{z a/ 



where z G C and we denote D := 4-. By substituting (BO) into 



obtain 



2/3 — Z a 

d_ 

df 



and 



1 c -D(S(z)-S(z a )) _ 1 e -Q(5( Zl )-g( Za )) 

+ _J_ e -^(*0(sw-s(*.)) = o, zee, z-l e G, 

Z\ — 2 a 



i e -D(S(*)-S(zfj)) _ 1 e -D(S(g2)-S(z p )) 



(81) 



+ _J_ e -^(*)(sW-s(*/»)) = , z g C, 2 2 G G. 



(82) 



2 2 - 2/3 



These identities are generating equations for (|69|). Indeed, expanding their 
left-hand sides in Taylor series in Z\ and 22, one gets (p9f). Furthermore, 
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evaluating flSTD at z = z@ and z = z a , one deduces that S(z a ) — S(zp) = 
DF a = DFp and so on. Hence 

F a = Fp = F = liiT d2 j mj . (83) 
Thus there is only one r-function for the d2DTL hierarchy. This feature is 



in agreement with earlier results ||11|| . Considering (pl|) at z = z\ E G and 
z = z 2 G G, and using flS2|), one shows that the function F satisfies 



(84) 



(z, - ^ - ^)e^^ - [z x - zp){z x - z a )e D ^ F + 
+ (z a - z p )(z! - Zl )e Da ^ Va ^ F = 0, 

and 

(z 2 - z a )(z 2 - z p )e- DV ^ F - (z 2 - z a )(~z 2 - z p )e- DV ^ F + 
+ {z p - z a )(z 2 - Z2 )e v ^ v ^)F = 



(85) 



These identities imply, in particular, the addition formulae for the dKP hi- 
erarchy. 

Then, evaluating fl8T|) at z — z 2 , using the relations S(z 2 ) = S(zp) — 
Vp(z 2 )F and S(z a ) — S(zp) = DF, one obtains 

1 + (Zg ~ Zp)(Zi - Za) c V a ( Zl )V fi ( Z2 )F 

{zi - z )(z 2 — z a ) 
_ {z\ - z a )(z 2 - zp) ^ Vg(z2) __ Va{zi)+D]DF _ ^ 
{zi -zp)(z 2 - z a ) 



Evaluating equation ( [82] ) at z — Z\ and taking into account that S(zi) = 
S(z a ) — V a (zi)F , one gets the same equation fl86l) . 

Equations (|S4"D, ( |85| ) and (|86| ) form a complete set of addition formulae 
for the d2DTL hierarchy which completely define [11], [18] its associated r 
function. 

Similarly to the one-point case, one can derive a generating equation for 
the hierarchy of equations for S^see p7| ) 

Finally, let us include into consideration the variations of positions z a and 
zp of singularities. The simplest case is given by the one logarithmic source 

5 = t ln^. (87) 

Z-Z/3 
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Since 

$z a = h ~ — , S z = h - — , (88) 

Z — Z a OZ a Z — Zp OZp 

one obtains the following equations (we use the gauge S(zo) = 0) 

S Za = Ue~ St - U, S Zp = Ve s * - V, (89) 

where 

U = e §tiZa \ V = — e~^\ (90) 

Equations (|89|) imply that 9 = lg(UV) obeys 

6 ZaZf) + (e% = 0, (91) 



which is again the d2DTL equation ([731) , but now the spacial variables are the 
positions of the logarithmic singularities. In particular for z a — w, zp = w 
one gets the elliptic version of the d2DTL hierarchy. Despite of the fact that 
(|9l| ) coincides with (|73|) their corresponding dressing procedures are different. 



Indeed, for the simplest choice S(z a ) := ip = 0, S(zp) := = 0, the formula 
flT5| ) gives the trivial solution 6 = const, of (|75[), while for (^TJ) one gets the 



nontrivial solution 

t 

.Zp — z a - 

For the whole d2DTL hierarchy one has two relations of the type 



e = \g(uv) = 2\ g (—!—). 

\Zr — Z n ' 



5 Multi-point sources 

For the general N points-like source 

N N oo 



a=l a=l n=l ^ " 



-a j 



one has N hierarchies of equations 

n 

St na = J2 U ^ e ~ kSt ° a > a = l,...,N,n>l. (93) 



fc=0 
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These families of equations are related by equations of the type (|5"E| ) 

e -St Oa s to0 + Aa(3e St 0a + ^ Q e" 5 ^ + 5 Q/3 = 0, (94) 

where 

A aP = e- St ^ iz «\ B a/3 = -(1 + A a/3 + Ap a ). (95) 

The hierarchy of equations (0), (@) is nothing but the universal Whitham 
hierarchy introduced in ||. It contains a number of interesting subhierar- 
chies. Obviously, there are N d(KP-mKP) hierarchies associated with each 
of the singularity points z a (a = 1, . . . , N). Then there are d2DTL hierar- 
chies corresponding to each pair (z a ,Zp) of the singularity points. One can 
show that there is only one function F for ([93| ) , (|9~4]) || .To demonstrate this 
property it is enough to apply the same argumentation as that used in the 
previous section to identify each pair F a and Fp. 

We will discuss now some of the subhierarchies of (|93|) ,(|94] ) with more 
that two singularity points. Firstly, by restricting ourselves to the variables 
tx a (a — 1, . . . , N), we obtain the following set of equations (in the gauge 
S(z ) = 0) 

%^-^%-^l%=0, a^P,a,P = l,...,N, (96) 

where £ Q := t\ a and ip a := S(z a ). By evaluating the left-hand side of fl96|) at 
the points z = z 7 (7 ^ ®,/3) one gets the bilinear system of equations 

vlpl, - - 4yl = 07) 

where a, (3 and 7 are different. Equivalent ly 

.a 

4^ + ^ = 1. (98) 

Equations ([M]) and ([)?]) can be treated as the dispersionless limits of the 

systems 

= {\gH a ) u ^ Ua + (lgtf^Vv (99) 
Hl U0 = (kH a ) Up Hl + (\gH% a H^ (100) 
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where a, (3 and 7 are different. The equations ( |100| ) form the Darboux system 
describing the iV-conjugate systems of surfaces in Euclidean space, and ( |99| ) 
are the corresponding equations for the position vector ||42|| . In this case the 
slow variables are £ a = eu a and the quasiclassical limit is implemented by the 
expressions ip = exp -, H a = exp — . Thus, one shows that (|9^) anbd (|100| ) 
represent the dispersionless limit of the Darboux system and its associated 
linear system (see also f37j). 

Our second example is determined by an So function of the form 

N M , 

s = ^—+yY,^+ t Y,-3^> ( 101 ) 

z z n=i z z n m=i z z m 

where a n , b m are arbitrary complex parameters and the sets {z n } and {z m } 
are disjoint. In this case 

P := S x = (- S x , 

z- z 

N 

a r 



(102) 



J yi 

. — ... 

n=l 

z z rn 

m=l 



From these expressions the following Hamilton- Jacobi type equations follow 

N V M V 

Sy = V + J2-^-, S^Vo + j:-^, (103) 

n=1 Ox Un m=1 b x - U m 

where 

U n = S x (z n ,x,y,t), U m = S x (z m ,x,y,t), (104) 

and Vo, V n , Vo, V m are functions of (x,y,t). The type of equations ( |103|) as 
well as its associated integrable systems have been considered in [ |T0|j . 
Another interesting subhierarchy corresponds to the choice 

, N 

So = -^ + T ^— 2 +yY,^^- ( 105 ) 

Z-Zq {z-Zq) 2 ^[Z-Z n 



25 



In this case one has the following equations 

- V 

Sy = Vo + J2lHyr> S t = S 2 x + U , (106) 

n=l 

where U n := S x (z n ,x,y,t). The associated integrable system 

u n , t + U 0>x - {U 2 n ) x = 0, 

V n>t - 2{V n U n ) x = 0, 1 < n < N, 

n (107) 

U , y + 2j2Vn,x = 0, 

i=l 

is the generalization of the Benney system to 2 + 1 dimensions proposed in 
1,0 

Other particular cases of the hierarchy (0),(@) will be considered else- 
where. 

Remark 5.1 All of our constructions in this and previous sections were 
local ones, with singularities (sources) located in certain domains of C. This 
means that one will get the same formulae by considering instead a chart of 
the Riemann sphere. Compactness of the Riemann sphere imposes certain 
constraints of the form ([92]) of the sources (singularities). In particular, 
E«=i*oa = (seee.g.H). 
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